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Abstract 

One of the fundamental properties of the quark-gluon plasma (QGP), the equation of state, is a 
subject of extensive studies in lattice QCD and an essential requirement for the correct hydro- 
dynamic modeling of heavy-ion collisions. Lattice QCD provides first-principle calculations for 
the physics in the non-perturbative regime. In this contribution, we report on recent progress by 
the HotQCD collaboration in studying the 2+1 flavor equation of state on lattices with the tem- 
poral extent N T = 6, 8, 10 and 12 in Highly Improved Staggered Quarks (HISQ) discretization 
scheme. Comparisons with equation of state calculations with different fermion actions are also 
discussed. 

At fixed cutoff (in our case, lattice spacing, a) the trace of the energy-momentum tensor, or 
interaction measure, can be related to the partition function of the system: 
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where s is the energy density and p is pressure. 

The equation of state, i.e., the dependence of pressure on temperature can be expressed then 
as an integral: 
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To calculate the interaction measure in 2+1 flavor QCD we work on a Euclidean space-time 
lattice and evaluate path integrals with the importance sampling technique. We use a tree-level 
improved action for gauge fields and the Highly Improved Staggered Quarks (HISQ) action for 
fermions [1]. To get finite values, subtraction of UV divergences is required: 
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where the subscript "0" refers to T = and "T" to finite temperature for observables evaluated 
at the same value of the cutoff. On a hyper-cubic lattice xA^ T the physical temperature is set by 



'The HotQCD collaboration members are: A. Bazavov, T. Bhattacharya, M. Buchoff, M. Cheng, N. Christ, C. DeTar, 
H.-T. Ding, S. Gottlieb, R. Gupta, P. Hegde, U. Heller, C. Jung, F. Karsch, E. Laermann, L. Levkova, Z. Lin, R. Mawhin- 
ney, S. Mukherjee, P. Petreczky, D. Renfrew, C. Schmidt, C. Schroeder, W. Soeldner, R. Soltz, R. Sugar, D. Toussaint, 
P. Vranas. 

Preprint submitted to Nuclear Physics A October 24, 2012 




Figure 1 : The difference in temperature (left) set by r\ and fa scale, described in the text. The running of the gauge 
coupling (right). 



the size of the temporal dimension and the lattice spacing as T — l/(N T a). For T — calculations 
we use N T ^ N s , and for T > we keep N S /N T = 4 and at fixed N T vary the lattice spacing a 
by varying the gauge coupling £1 = 10/g 2 . The continuum limit in this setup is controlled by 
N T — » co, therefore, we carried out this study on several sets of ensembles with N T = 6, 8, 10 
and 12. The strange quark mass m s is tuned to the physical value, while the two degenerate 
light quarks have masses mi = m s /2Q, slightly heavier than physical (mi - m,/27). To set the 
lattice spacing in physical units (fm) we use the Sommer-type scale |2] n = 0.3106 fm or, 
alternatively, the kaon decay constant, = 156.1 MeV. From Fig. Q] (left) one can see what 
effect using different reference observables has on converting the temperature from lattice units 
to MeV. Over the temperature range of interest on N T = 6 lattices the difference is within 9 MeV, 
and on N T = 12 within 2 MeV. Another effect of using different scales is related to running of the 
corresponding /^-functions: 

dB 1 dm s (B) 

Rfi(P) = -a-f, R m {p)= — -gp* , (4) 
da ni s (p) dp 

where m s (J3) defines a line of constant physics (LCP), i.e., such combination of the gauge cou- 
pling and the strange quark mass so that pion and kaon masses (in MeV) stay approximately 
constant in the whole (3 range used in the simulation. Rp for r\ and scale is shown in Fig.Q] 
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Figure 2: Comparison of the HISQ/tree interaction measure on N T = 6, 8, 10 and 12 lattices with the stout continuum 
estimate |4] (left), N T = 6 and 10 HISQ/tree data (right). Filled (open) symbols in the right panel correspond to the r\ 
(Jk) scale, see text. 



(right), together with the 2-loop perturbative result. As is clear from Eq. (O, the /3-functions 
enter into the interaction measure multiplicatively. 

Present status of the interaction measure with HISQ/tree is shown in Fig. [2]together with the 
stout continuum estimate of Ref (4J] . Although we have not yet performed the continuum limit, 
the HISQ/tree results seem to disagree with stout in T — 170 - 350 MeV range. 

The interaction measure for N T - 6 and 10 is shown in Fig.[2](right), N T — 8 and N T = 12 in 
Fig. [3] We also include our data with the asqtad action, where available. For comparison we plot 
the HISQ/tree and asqtad results with the temperature scale set with both r\ (filled symbols) and 
fx (open symbols). 
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Figure 3: The interaction measure for the HISQ/tree action on N T = 8 (top) and N T = 12 (bottom) ensembles with r\ 
(left) and (right) scale. 

On the coarsest lattice, A^ T = 6, Fig. [2] (right) the cutoff effects for the HISQ/tree action are 
easily observed: the data plotted with fx scale (red open symbols) lies lower and is shifted to the 
left, compared to the data with the r\ scale (red filled symbols). On N T = 8 the cutoff effects for 
HISQ/tree are milder, but still are very substantial for the asqtad action, compare the filled blue 
symbols in the left panel to the open blue symbols in the right panel of the upper part of Fig. [3] 
On N T = 10 and 12 for both HISQ/tree and asqtad, Fig.|2](right) and the lower part of Fig. |3] the 
cutoff effects are comparable with the statistical errors and agreement between the two actions is 
good for r\ and fx scale. 

The low-temperature behavior of the interaction measure is shown in Fig. |4] (left). The lines 
represent the hadron resonance gas (HRG) calculation with the physical pion and also with the 
full unphysical heavier pion multiplet, as encountered on the lattice. (The effects of heavier 
pion multiplet at fixed lattice spacing for the HISQ/tree action are quantitatively addressed in 
Ref. There are two interesting observations: 1. Although the root-mean-squared pion mass 
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varies from 400 to 200 MeV for N T - 6 - 12 ensembles, the interaction measure in HRG is not 
very sensitive to this variation. Thus, better agreement of the HISQ/tree lattice data for various 
N T in the low-temperature region is due to both, improved properties of the HISQ/tree action and 
little sensitivity of this observable to the cutoff effects in the pion sector. 2. The lattice result 
starts to disagree with the HRG model at T ~ 150 - 160 MeV. 

The interaction measure at high temperature is presented in Fig. [4] (right) for the HISQ/tree, 
p4 and stout action. The lines represent perturbative calculations. For T up to 800 MeV cutoff 
effects expected in the T — > oo limit are observed only for the p4 action. 



1. Conclusion 

We continued the calculation of the 2+1 flavor QCD equation of state with the HISQ/tree 
action on N T = 6, 8, 10 and 12 ensembles. Using different observables to set the scale, r\ and 
/k, allows for a crude estimate of the magnitude of cutoff effects, which seem to be small at the 
finest N T = 12 lattices. Better control over the statistics on N T = 10 and 12 ensembles is required 
before an extrapolation to the continuum may be attempted. 
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